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Abstract 

Companion stars exposed to high energy radiation from a primary neutron star or 
accreting black hole can experience significant spallation of their heavy elements, so that 
their atmospheres would be extremely rich in lithium, beryllium, and especially boron. 
In this paper we note that the detection or non-detection of these elements, and their 
relative abundances if detected, would provide a diagnostic of the high energy output of 
the primary, and possibly the shock acceleration of particles at the companion's bow shock 
in a pulsar wind. 
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1. Introduction 

Pulsar emission theory is by now an established if unresolved topic. Radio pulsations 
are believed to imply counterstreaming pairs, which are produced via the curvature ra- 
diation of very high energy gamma rays. Below the "death line", it is believed that the 
primary radiation is mainly curvature radiation that fails to develop a next generation 
of pairs. Even while pulsing, a pulsar could easily put out ~ 10~^ of its spin down en- 
ergy into pulsed gamma rays, (Usov 1983) since a) pair production shorts out the field to 
the point that the pairs are accelerated only enough to marginally produce pairs, and b) 



charged particles that are accelerated outward will eventually reach a site where emitted 
photons can escape without further pair production. A reasonable estimate for curvature 
radiation is probably the product of the Goldreich-Julian current and the polar voltage 
drop, whence the fraction of total spin down power radiated as curvature gamma rays can 
be of order 10~^ or more for millisecond pulsars. Luminous, nearby pulsars such as the 
Crab are observed to yield pulsed gamma rays at a total power consistent with theoretical 
expectations. Recent EGRET limits are also consistent with these expectations (Fierro et. 
al. 1995). However, additional gamma radiation may be generated by pairs striking and 
cascading in the atmospheres of companions to pulsars. Also, the pair luminosity itself 
may be enhanced by shock acceleration of pulsar wind particles at the bow shock of its 
companion ( Arons and Tavani, 1994, Grove et. al. 1995). 

In this paper, we suggest that the output in high energy (E > 20 Mev) quanta from 
pulsars with close binary companions can be diagnosed via the production of light ele- 
ments in the companion's atmosphere by photospallation. The idea of using light elements 
as a dosimetric diagnostic of pulsars is not new: it has been discussed in the context of 
very young pulsars irradiating early supernova remnants (Eichler and Letaw 1987). More 
recently, Li has been reported to be present in roughly cosmic abundances in V404 Cygni, 
contrary to expectations that it should have been mostly destroyed. It has been considered 
that the Li has been replenished by photospallation caused by irradiation from its com- 
pact (primary) companion. Although this invokes some numerical coincidence, alternative 
explanations for the Li abundance in V404 Cygni-like systems are briefly mentioned below. 

The question of whether "black widow" type evaporation scenarios occur for accreting 
neutron stars and pulsars is still an open one. It has been suggested that heat induced 
evaporation can be important for companions to both accreting neutron stars and pulsars 
(Eichler and Ko 1988; Ruderman et aZ.1989; Ruderman, Shaham and Tavani 1989) though 
the efficiency of mass loss and extent of ablation have been questioned on both obser- 
vational and theoretical grounds (Levinson and Eichler 1991; Eichler 1991; Gedalin and 
Eichler 1993). This matter is particularly questionable for low luminosity pulsars where 
some mechanism must be invoked for converting the spin down power to a form of energy 
suitable for mass evaporation. In the latter case, some scenarios assume that the pulsar 
spin down power is somehow converted to soft gamma rays with high efficiency (Kluzniak 
et aZ.1988; Phinney et a/. 1988). Several eclipsing pulsars have been discovered at this time 
(Fruchter et aZ.1988; Lyne et aZ.1990; Johnston et a/. 1992) and the eclipse is suggestive 
of some mass loss mechanism, though not necessarily implicative of significant ablation. 
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In the PSR 1957+20 system, for example, viable eclipse mechanisms typically require a 
plasma frequency at the eclipse site of about 0.1 of the frequency of the radio waves that are 
being eclipsed (Gedalin and Eichler 1993). Assuming that the outflow proceeds at about 
the orbital velocity or somewhat higher, this implies a mass loss rate of about lO^^g/s. 

2. Photospallation and light elements 

The photospallation cross section for C^^ +7 B^^ +p becomes significant at photon 
energies above 15 Mev or so. It averages about 5 millibarns in the photon energy range 20 
to 25 MeV, or about 1 millibarn per unit logarithm in photon energy at 20 MeV (Taran 
and Gorbunov 1967). Higher energy gamma rays or pairs hitting the companion surface 
cascade via bremsstrahlung pair production cycles, and always pass through this energy 
range. Shock acceleration of pairs followed by synchrotron emission would convert much 
of the shock energy to photons of about several Mev, so the amount of photospallation 
predicted by this scenario depends sensitively on a detailed calculation of the synchrotron 
spectrum. 

We consider the consequences of the hypothesis that some fraction e of the pulsar's spin 
down power L arrives in the form of (or is converted to) 20 Mev photons. For convenience, 
assume the average photon energy included in this fraction is 10~^ erg. Then the hard 
photon flux incident on the companion is / = 3 x 10^'*/i4.5, where /14.5 = ( 3x10-3 
photons/cm^ s, L35 is L in units of 10^^ erg s~-^, and Dn is the orbital separation in units 
of lO-^-*^ cm. The lifetime of a carbon nucleus exposed to f is about 3 x 10-^^//i4.5 seconds. 

Can the top radiation length of material remain on the surface for more than 10^° 
seconds (in which case the heavy elements would be completely destroyed)? The mass loss 
rate for the PSR 1957+20 system is conservatively estimated to be of order 10^^^/s, and 
could easily be lO^^g/s (so we define it to be M13 times this amount), from a surface area 
of 167ri?2 X 10^°cm^, (here R2 is the radius in units of 2 x 10-^°cm. R2 = 1 corresponds 
to the companion to PSR 1957+20 filling its Roche lobe ) or a stripping rate of at least 
2 X 10~^Mi3i2^^ g/cm^s. As a radiation length is roughly 10^ gm/cm^, the stripping 
time could in principle exceed the lifetime of heavy nuclei at the surface, if there is no 
significant mixing of the surface layers, but it does not appear to exceed it by such a large 
factor that the spallation products themselves would be entirely broken up to helium and 
hydrogen. Moreover, the observed absorption lines, insofar as they affect the continuum, 
are consistent with standard solar abundances (Romani, private communication). We thus 
conclude that for PSR 1957+20, and similar system, the spallation of heavy elements, if it 
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occurs, does not defeat itself by the total destruction of either the heavy elements or the 
products. [Such destruction could in principle be important for some hypothetical range of 
parameters, (high /14.5, low M, requiring suppressed mass evaporation by the radiation). 
In this case however the destruction of primary heavy elements such as carbon would be 
the more conspicuous effect.] 

Let us first suppose that all of the spallation products that are produced are eventually 
evaporated. This implies that the exposure time at the surface is set by the mass loss rate, 
independent of convective mixing. By the above, the fraction of carbon that is spalled to 
lighter nuclei can be of order 1.5xlO~^/i4.5(M/10^^(7m/s) ^i?!? unity, whichever is less. 
For PSR 1957+20-type parameters, this is many orders of magnitude above the cosmic 
values of Li, Be and B. The largest change is in the abundance of boron. The abundance 
levels depend on the cross-sections of the reactions that produce and destroy the respective 
nuclei, the photon flux and the exposure time. In the case of spallation that is just below 
the limit of total destruction of heavy nuclei, the ratios B^'^+^^/H, Be^/H and Li^/H could 
be as high as 5.8 x 10""^, 5.9 x 10~^, 3.7x 10~^ (by number) after spallation (Boyd and Fencl 
1991) compared to the cosmic ratios of 3.4 x 10"^°, 4.5 x 10-^^ 2.3 x 10"^ respectively 
(Cameron 1982). We have estimated these ratios for a photon number spectrum with an 
index —2.7, assuming that the expected spectrum of 1957+20 will be similar to that of 
Crab nebula (Nolan et a/. 1993). These ratios change only slightly with the exact shape 
of the spectrum above 10 MeV. In any case, the spallation products should be above 
the otherwise expected levels, if the incident flux of sufficiently hard gamma radiation 
on the companion is a small fraction (10~^ or more for 1957+20) of the intercepted spin 
down power. We suggest HST observations of the BI resonance lines near 2497A which 
have been used to estimate the boron abundances in halo stars (Duncan, Lambert and 
Lembke 1992). Unfortunately, the reddenning of 1957+20 may make observations in the 
UV difficult, although this does not rule out the UV observations of other sources. The Li 
I resonance line at 6707. 8 A can be used to measure the lithium abundance (as was done 
for V404 Cygni, see below). 

The above neglects both dilution due to turbulent convection to deeper levels, which 
may occur more rapidly than the stripping time, and destruction by subvection to the 
core. In discussing this, it is useful to note that the column density of the entire star is 
M/7ri?^, where M and R are respectively the mass and radius of the star, is about 10^ 
radiation lengths. // the mass evaporation is such that all of the layer down to which the 
mixing obtains eventually is lost, then in steady state, the mixing has no effect on the level 



4 



of spallation products at the surface. This includes the case in which the entire star is 
eventually evaporated. If, on the other hand, the mixing depth is larger than the time 
integrated ablation depth A^ (which grows linearly in time), then the dilution of spallation 
products is simply the ratio of the two A^/A^. While an evaluation of the mixing depth is 
beyond the scope of this paper, we note that for PSR 1957+20, the layer of mass loss can 
be conservatively estimated from the dispersion measure of the eclipsing wind as lO^^^r/s 
times the age t. For t of order a spin down time, which for this particular system is roughly 
a Hubble time or more, the lost mass is then about 10~^ of the companion mass, so that 
the astrophysical uncertainty imposed by our ignorance of turbulent mixing is not all that 
large. Even if the system is only ~ 3 x 10^ yr old, the uncertainty in Xm/^a is several 
orders of magnitude less than the maximum light element enhancement. 

The light element abundance in the gas that has been blown off from the companion, 
is simply the ratio of the rate of accumulation of light element nuclei to the stripping rate, 
apart from the dilution factor due to convection as discussed in the previous paragraph. 
The rate of accumulation of light elements due to spallation has been calculated by Boyd 
and Fencl (1991). In the limit that neither carbon nor lithium destruction is significant at 
the surface we can use their calculations to express the abundances of light elements (for 
the case of A^ > A^) as: 

|^(-^)io-«Mr/i?i/M.5, 

^~(^)9xlO-^3Mr3^i?^A4.5, (1) 

where, as argued in the preceding paragraph, the ablation depth A^ is a linear function 
of time and depends on the particular pulsar and its companion. Eqns(l) show that the 
resulting B abundance in the companion can be very important even in the case of deep 
mixing, much more than those of Li and Be. We note here that the ratio Aa/A^ for a 
companion mass of 0.02 M0 can be written as O.lrioMisi?"^, where Rm is the mixing 
depth in units of i?2, and tiq is the system age in 10-^°yr. 

We have considered the possibility that spallation products are destroyed by nuclear 
burning. Probably the companion star is too cool for this, if it fills its Roche lobe, in the 
case of 1957+20, where the companion mass is only 0.025 and the radius about 0.3 (in 
solar units) (Fruchter et aZ.1988). The central temperature is for these numbers only about 
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1.5 X lO^K. In general, a light companion still on the main sequence could burn its spallation 
products; a star of 0.08 solar masses, for example, has a central temperature of 4 x lO^K, 
enough for lithium burning, which requires only 3 x lO^K (Swenson, Stringfellow and 
Faulkner 1990), but, clearly even a modest amount of bloating would cool the companion 
core too much to incinerate the spallation products. The case against incineration is even 
stronger for Be and B, and/or lighter companions, such as in the case of 1957+20, where 
the central temperature may not in general be enough for incineration under any set of 
assumptions. 

3. V404 Cygni 

The lithium abundance in the secondary in the V404 Cygni system has been estimated 
to be of the order of Li/H r-^ 10"^, which is close to the cosmic interstellar ratio (Martin 
et al. 1992). This is anomalous considering the fact that the secondary is most probably a 
G9 dwarf star (Martin et a/. 1992), and a G/K main sequence star should deplete its lithium 
abundance by an order of magnitude after ~ 100 Myr and by 2-3 orders of magnitude 
after ~ 1 Gyr. In the light of the above discussion on photospallation, it seems possible 
that 7-radiation from the vicinity of the central object, a stellar mass black hole, can 
easily reimburse the lost lithium in the secondary star, as has been noted by previous 
authors (Martin et a/. 1992, Gilmore 1992; also see discussions on photospallation and light 
elements in other astrophysical contexts, e.g., Boyd and Fencl 1991, Gnedin and Ostriker 
1992) (While Martin et aZ.(1992) proposed particle spallation in the accretion disk around 
the primary, spallation at the companion's surface should work just as well). The observed 
abundance of lithium that is near the cosmic abundance would require a rough, coincidental 
balance between the production and destruction processes (at the present epoch) described 
in the eqns (1). However, there may be two other reasons for the preservation of Li at 
roughly cosmic abundances that do not require such a coincidence: Firstly, it may be that 
heating of the companion surface by the primary compact object merely stabilizes the outer 
layers to convection, so that the light elements in these layers are never subducted to the 
depths at which they are destroyed. Secondly, it may be that the black widow effect, if the 
companion is destroyed on a timescale of ~ 10 Myr, selects out for observation only those 
systems that are too young to have suffered significant Li destruction. Theoretical analysis 
of the possibilities from first principles, in our view, are extremely difficult, and not in any 
case the primary topic of this paper. Their viability depends, as for the case of pulsar 
companions, on the turbulence induced by horizontal pressure gradients resulting from 
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uneven heating of the surface. Such turbulence, together with heat conduction, competes 
with outward mass flow in determining the entropy distribution in the outer layers of 
the heated companion, and the matter is in our view not fully resolved. But if the two 
"preservation" scenarios for the light elements can be observationally distinguished from 
the scenario of "restoration" by photospallation (the former predicting a higher ratio of Li^ 
to the other light elements), then the results might have implications for the companions 
to pulsars as well. 

Conclusions 

Any model of companion irradiation by a primary pulsar in which quanta exceeding 
15 Mev or so strike the companion surface predicts that the companion atmosphere be 
anomalously rich in Li, Be and, particularly B if the central temperature of the companion 
is too small to burn it. Spectroscopic observations of an ever growing number of compan- 
ions to binary pulsars may prove to be a powerful diagnostic of the high energy output of 
the primary pulsar that is predicted theoretically. 

We are grateful to Profs. V.V. Usov and R. Romani for useful discussions. This 
research was supported in part by the Israeli Foundation for Basic Research. 
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